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ABSTRACT 

The  transport  and  mixing  characteristics  of  a  large  sample  of  air  parcels  within  a  mature  and  vertically 
sheared  hurricane  vortex  are  examined.  Data  from  a  high-resolution  (2-km  horizontal  grid  spacing)  nu¬ 
merical  simulation  of  real-case  Hurricane  Bonnie  (1998)  are  used  to  calculate  Lagrangian  trajectories  of  air 
parcels  in  various  subdomains  of  the  hurricane  (namely,  the  eye,  eyewall,  and  near  environment)  to  study 
the  degree  of  interaction  (transport  and  mixing)  between  these  subdomains.  It  is  found  that  1)  there  is 
transport  and  mixing  from  the  low-level  eye  to  the  eyewall  that  carries  air  possessing  relatively  high  values 
of  equivalent  potential  temperature  (0e),  which  can  enhance  the  efficiency  of  the  hurricane  heat  engine;  2) 
a  portion  of  the  low-level  inflow  of  the  hurricane  bypasses  the  eyewall  to  enter  the  eye,  and  this  air  both 
replaces  the  mass  of  the  low-level  eye  and  lingers  for  a  sufficient  time  (order  1  h)  to  acquire  enhanced 
entropy  characteristics  through  interaction  with  the  ocean  beneath  the  eye;  3)  air  in  the  mid-  to  upper-level 
eye  is  exchanged  with  the  eyewall  such  that  more  than  half  the  air  of  the  eye  is  exchanged  in  5  h  in  this  case 
of  a  sheared  hurricane;  and  4)  that  one -fifth  of  the  mass  in  the  eyewall  at  a  height  of  5  km  has  an  origin  in 
the  mid-  to  upper-level  environment  where  6e  is  much  less  than  in  the  eyewall,  which  ventilates  the 
ensemble  average  eyewall  Be  by  about  1  K.  Implications  of  these  findings  for  the  problem  of  hurricane 
intensity  forecasting  are  briefly  discussed. 


L  Introduction 

Vortices  exhibit  the  well-known  property  that  mixing 
between  their  cores  and  near  environment  is  strongly 
suppressed.  Examples  are  manifold,  including  vortex 
rings  (Vladimirov  and  Tarasov  1979),  mesoscale  ocean 
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vortices  (Robinson  1982),  emergent  vortices  in  quasi¬ 
geostrop  hie  and  two-dimensional  turbulence 
(Carnevale  et  al.  1991;  McWilliams  et  al.  1994),  the 
wintertime  stratospheric  vortex  (McIntyre  1989),  and 
tropical  cyclones  (TCs;  Willoughby  1998).  The  details 
of  these  flows  are  quite  different,  yet  the  tendency  for 
vortices  to  act  as  mixing  barriers  to  environmental  fluid 
is  believed  central  to  both  their  emergence  and  long 
lifetimes  in  complex  flows  (e.g.,  McWilliams  1984;  Mi- 
zuta  and  Yoden  2001).  This  vortex  robustness  can  be 
traced  to  the  quasi-elastic  behavior  of  vortices,  which 
results  from  having  mechanisms  that  can  oppose  distur¬ 
bances  that  are  either  axisymmetric  [via  centrifugal  res¬ 
toration  (Howard  and  Gupta  1962;  Charney  1973)]  or 
nonaxisymmetric  [via  realness  of  discrete  spectra 
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(Vladimirov  and  Tarasov  1980)  or  via  a  tendency  to 
axisymmetrize  (Melander  et  al.  1987;  Montgomery  and 
Kallenbach  1997;  Bassom  and  Gilbert  1999;  Reasor  et 
al.  2004;  Schecter  and  Montgomery  2004)].  If  the  ro¬ 
bustness  of  a  vortex  core  were  to  be  reduced,  however, 
this  could  have  important  consequences  on  the  main¬ 
tenance  of  that  vortex. 

In  TCs,  the  manifestation  of  some  compromise  in  the 
robustness  of  the  vortex  is  the  stirring1  of  air  (and  at¬ 
tendant  properties  such  as  angular  momentum  or 
equivalent  potential  temperature  6e)  between  the  vari¬ 
ous  domains  of  the  storm;  that  is,  the  eye,  the  eyewall, 
and  the  surrounding  environment.  Among  the  key  at¬ 
mospheric  mechanisms  that  have  been  identified  to 
produce  mixing  between  TC  domains  are  transport  by 
frictionally  induced  inflow  that  slips  under  the  eyewall 
(Persing  and  Montgomery  2003,  hereafter  PM03;  Braun 
2002);  stirring  by  barotropic/baroclinic  (inflectional)  in¬ 
stabilities  that  break  down  the  high  potential  vorticity 
annulus  comprising  the  eyewall  (Schubert  et  al.  1999; 
Kossin  and  Schubert  2001;  Nolan  and  Montgomery 
2002;  Kwon  and  Frank  2005),  and  stirring  by  an  im¬ 
posed  vertical  wind  shear  (Frank  and  Ritchie  2001; 
Zehr  2003;  Emanuel  et  al.  2004).  Each  of  these  pro¬ 
cesses  has  been  shown  to  be  capable  of  significantly 
altering  the  intensity  of  TCs  (Wu  and  Braun  2004;  Kos¬ 
sin  and  Eastin  2001;  Bender  1997;  Wang  2002;  Emanuel 
1989,  1997;  Montgomery  et  al.  2006;  Frank  and  Ritchie 
2001;  Wong  and  Chan  2004).  While  signifying  a  com¬ 
promise  in  the  robustness  of  a  vortex,  these  mixing  pro¬ 
cesses  can  alter  the  thermodynamics  of  the  hurricane 
eyewall,  whether  by  enhancement  or  dilution  of  rel¬ 
evant  thermodynamic  properties.  Recall  for  example, 
that  the  uptake  of  moist  enthalpy  from  the  underlying 
ocean  and  the  ascent  in  the  eyewall  of  the  TC  represent 
the  isothermal  and  moist  adiabatic  legs  of  an  idealized 
Carnot  cycle,  respectively.  An  enhanced  intake  of  en¬ 
ergy  to  the  Carnot  cycle  at  this  stage  suggests  an  in¬ 
crease  in  useful  work  (enhancement  of  the  kinetic  en¬ 
ergy  field)  and  ultimately  an  increase  in  the  loss  of 
energy  to  frictional  dissipation  at  the  ocean  surface 
(Emanuel  1986,  1995;  PM03). 

Evidence  in  support  of  the  hypothesis  that  instabili¬ 
ties  in  the  TC  eyewall,  and  their  finite  amplitude  co- 


1  We  follow  the  definitions  found  in  Haynes  (2002).  “Trans¬ 
port”  is  to  carry  air  with  various  properties  from  their  source  to 
their  sinks.  “Stirring”  is  to  bring  air  with  different  properties  into 
closer  proximity  with  each  other.  Stirring  that  proceeds  to  the 
point  that  the  originating  characteristics  of  the  air  can  no  longer 
be  distinguished  is  called  “mixing”.  Mixing  will  also  be  used  in 
general  discussion  where  the  distinction  between  these  processes 
is  not  important. 


herent  structures  [such  as  illustrated  in  Schubert  et  al. 
(1999)  and  Montgomery  et  al.  (2002)],  contribute  posi¬ 
tively  to  the  intensity  of  the  storm  was  presented  by 
Emanuel  (1989,  1997)  and  PM03.  In  Emanuel  (1997), 
the  breakdown  of  the  vortex  sheet  comprising  the  eye- 
wall  of  a  developing  storm  was  shown  to  be  essential  for 
achieving  the  maximum  potential  intensity  predicted  by 
the  Carnot  theory  (hereafter  E-PI).  In  PM03,  eyewall 
mesovortices  bordering  the  low-level  eye  and  eyewall 
region  were  spawned  by  a  (Kelvin-Helmholtz)  shear 
instability  on  the  annulus  of  high  toroidal  vorticity  as¬ 
sociated  with  the  secondary  circulation  of  the  hurri¬ 
cane.  A  characteristic  of  these  coherent  structures,  not 
anticipated  by  Emanuel  (1989,  1997),  was  their  ability 
to  persistently  flux  high  entropy  air  from  the  low-level 
eye  to  the  eyewall.  The  temporal  eddy  and  mean  trans¬ 
port  of  high  entropy  air  from  the  low-level  eye  to  the 
eyewall  was  shown  to  provide  supplementary  fuel  to 
the  heat  engine,  in  addition  to  what  is  drawn  from  the 
ocean  underneath  and  outside  the  eyewall.  By  utilizing 
this  extra  fuel,  the  hurricane  can  obtain  a  superintense 
state  (PM03;  Montgomery  et  al.  2006),  in  which  the 
storm  intensity  exceeds  E-PI.  Even  if  a  storm  in  nature 
is  not  superintense,  PM03  suggested  that  the  super¬ 
intensity  mechanism  may  still  operate  to  partially  miti¬ 
gate  the  negative  influences  associated  with  upwelling 
cooler  water,  cold  downdrafts,  environmental  vertical 
shear,  or  interactions  with  midtropospheric  dry  air. 

A  complementary  situation  to  superintensity  is  the 
problem  of  transport  and  stirring  of  midtropospheric 
dry,  low-entropy  air  from  outside  the  rain  area  into  the 
eyewall  region.  This  is  a  form  of  entrainment  that  oc¬ 
curs  on  the  vortex  scale.  Previous  work  (Simpson  and 
Riehl  1958;  Riehl  and  Malkus  1961;  Gray  1968;  Eman¬ 
uel  et  al.  2004)  has  suggested  that  such  an  exchange 
process,  which  can  result  from  an  imposed  vertical  wind 
shear,  can  dilute  the  heat  content  of  the  eyewall  thereby 
weakening  the  intensity  of  a  TC.  This  process  was 
coined  ventilation  by  Simpson  and  Riehl  (1958).  In 
practical  situations,  ventilation  may  occur  as  vertical 
wind  shear  impinges  on  a  storm  or,  alternatively,  as  a 
storm  approaches  a  mid-  and  upper-tropospheric  air 
mass  with  low  entropy.  During  the  process  of  producing 
an  intensity  forecast,  forecasters  must  decide  if  the  TC 
core  will  be  adversely  affected  by  this  interaction.2  In 
general,  ventilation  requires  dry  air  to  breach  the  in¬ 
tense  vorticity  of  the  eyewall.  Idealized  modeling  stud¬ 
ies  suggest,  however,  that  this  can  only  occur  via  non¬ 
linear  breakdown  of  the  annulus  of  high  potential  vor- 


2  In  this  paper,  the  core  of  the  TC  refers  to  the  eye  and  the 
high-vorticity  region  of  the  eyewall. 
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ticity  comprising  the  eyewall  and  the  attendant  stirring 
between  the  environment  and  eyewall  air  masses  (e.g., 
Schubert  et  al.  1999).  Ventilation  is  unlike  the  superin¬ 
tensity  mechanism  that  provides  extra  fuel  from  the 
low-level  eye  region.  Processes  that  entrain  dry  air  into 
the  TC  core  will  be  referred  to  here  as  the  anti-fuel 
mechanism. 

Little  is  known  quantitatively  about  the  efficiency  of 
the  fuel  or  anti-fuel  mechanisms  in  real-case  TCs.  In 
this  paper  we  take  the  first  step  in  a  more  systematic 
study  to  test  and  quantify  these  ideas.  The  real-case 
numerical  simulation  of  Hurricane  Bonnie  (1998)  by 
Braun  et  al.  (2006)  employing  the  nonhydrostatic  fifth- 
generation  Pennsylvania  State  University-National 
Center  for  Atmospheric  Research  (PSU-NCAR)  Me- 
soscale  Model  (MM5,  version  3.4;  Dudhia  1993;  Grell  et 
al.  1995)  provides  a  convenient  starting  point  for  our 
examination.  By  using  a  large  sample  of  air  parcels 
seeded  throughout  the  storm  we  build  upon  the  study 
of  Braun  (2002)  and  compute  several  thousand  three- 
dimensional  Lagrangian  trajectories  for  air  parcels 
seeded  throughout  the  eye,  eyewall,  and  surrounding 
regions  of  Hurricane  Bonnie.  The  large  number  of 
Lagrangian  trajectories  is  used  to  gather  a  census  of  the 
behavior  and  thermodynamic  properties  of  air  parcels 
in  various  regions  of  the  simulated  hurricane. 

Although  Bonnie  represents  only  a  single  case,  and 
results  gleaned  from  it  cannot  easily  be  generalized,  we 
believe  the  methodology  developed  herein  is  nonethe¬ 
less  useful  for  laying  the  groundwork  for  a  more  sys¬ 
tematic  study  of  the  hurricane  mixing  phenomena  using 
idealized  TC  flow  configurations.  Our  objective  here  is 
threefold:  Identify  the  nature  of  the  mixing  processes 
that  occur  within  this  vertically  sheared  storm;  quantify 
the  efficiency  of  the  mixing  between  the  eye,  eyewall, 
and  environment;  and  relate  the  current  findings  on 
mixing  within  the  vortex  to  the  heat  transport  that 
maintains  the  vortex  at  approximately  constant  inten¬ 
sity  (in  terms  of  maximum  tangential  wind)  during  a 
period  of  near-constant  vertical  shear. 

The  outline  of  the  paper  is  as  follows.  Section  2  de¬ 
scribes  the  model  setup  and  provides  an  overview  of  the 
Hurricane  Bonnie  simulation.  Section  3  describes  the 
Lagrangian  trajectory  seeding  strategy  and  provides  re¬ 
sults  illustrating  the  mixing  that  occurs  between  various 
regions  of  the  model  hurricane.  Figure  1  provides  a 
simplified  view  of  these  mixing  processes  in  the  context 
of  a  hurricane  that  experiences  a  moderate  degree  of 
vertical  wind  shear  [10  ms-1  (10  km)^1].  Section  3a 
focuses  on  the  low-level  eye/eyewall  interaction,  includ¬ 
ing  evidence  of  the  superintensity  mechanism  (trajec¬ 
tory  class  I  in  Fig.  1);  section  3b  focuses  on  trajectories 
in  the  boundary  layer  inflow  layer  and  a  quantification 


of  the  low-level  eye  mass  replacement  (class  II  in  Fig. 
1);  section  3c  focuses  on  the  upper-level  eye/eyewall 
interaction  and  investigates  the  eye  containment  hy¬ 
pothesis  proposed  by  Willoughby  (1998,  hereafter 
W98)  (class  III  in  Fig.  1);  and  section  3d  focuses  on  the 
interaction  of  the  midlevel  environment  with  the  storm 
and  how  the  environmental  vertical  shear  may  serve  to 
ventilate  the  eyewall  and/or  eyewall  region  (class  IV  in 
Fig.  1).  Section  4  concludes  with  a  summary  of  the  re¬ 
sults  and  objectives  for  future  research. 

2.  Model  setup  and  overview  of  the  Bonnie 
simulation 

Braun  et  al.  (2006)  presented  an  in-depth  analysis  of 
a  numerical  simulation  of  Hurricane  Bonnie  using  the 
PSU-NCAR  MM5  (version  3.4).  We  use  the  same 
simulation  in  this  study  but  focus  on  the  mixing  char¬ 
acteristics  that  can  be  deduced  from  the  use  of  a  large 
number  of  Lagrangian  trajectories. 

a.  Experimental  design 

Summarizing  the  numerical  setup,  four  nested  grid 
domains  of  increasing  horizontal  resolution  were  used. 
Because  of  computational  limits,  the  simulation  was 
carried  out  in  two  steps.  First,  the  two  coarser  meshes 
were  simulated,  using  a  36-km  horizontal  grid  spacing 
and  91  X  97  X  27  grid  points  in  the  x,  y,  and  z  directions 
for  the  first  grid,  and  160  X  160  X  27  grid  points  at 
12-km  horizontal  grid  spacing  on  the  second  grid.  The 
simulation  on  these  two  outer  grid  meshes  was  started 
at  1200  UTC  22  August  1998  and  run  for  36  h,  with 
model  output  saved  every  hour.  A  modified  version  of 
the  Blackadar  planetary  boundary  layer  scheme  is  em¬ 
ployed  in  which  surface  roughness  is  calculated  for  mo¬ 
mentum,  temperature,  and  moisture  following  Garratt 
(1992)  and  Pagowski  and  Moore  (2001).  Cloud  effects 
are  parameterized  on  the  coarse  mesh  in  two  ways:  the 
Grell  cumulus  scheme  (Grell  et  al.  1995)  is  used  to 
represent  unresolved  convective  processes  while  the 
Goddard  microphysics  scheme  (Tao  and  Simpson  1993; 
McCumber  et  al.  1991)  is  used  to  represent  grid-scale 
cloud  microphysical  processes. 

A  higher-resolution  simulation  was  then  carried  out 
using  the  1-h  output  from  the  36-  and  12-km  grid  results 
to  provide  initial  and  boundary  conditions.  The  model 
grids  consisted  of  a  6-km  grid  with  225  X  225  X  27  grid 
points  and  a  2-km  nest  with  226  X  226  X  27  grid  points. 
The  simulation  on  these  high-resolution  grids  was 
started  6  h  into  the  coarse-resolution  forecast  to  allow 
for  model  spinup  on  the  12-km  grid,  and  was  carried  out 
for  30  h  until  0000  UTC  24  August.  The  2-km  grid  was 
moved  periodically  to  keep  it  centered  over  the  storm. 
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Fig.  1.  A  simplified  schematic  of  a  moderately  sheared  hurricane.  The  dark  (light)  gray 
lines  show  the  approximate  outer  (inner)  edge  of  the  eyewall.  The  dotted  line  represents 
the  inversion  level  in  the  eye.  The  four  subsections  of  section  3  examine  different  mixing 
process  within  the  hurricane  using  different  trajectory  strategies  for  the  numerical  simu¬ 
lation  of  Hurricane  Bonnie  by  Braun  et  al.  (2006).  Section  3a  examines  the  mixing  from  the 
low-level  eye  to  the  eyewall  (trajectory  class  I)  by  counting  the  number  of  trajectories  that 
enter  the  eyewall  vs  those  that  remain  in  the  eye.  Section  3b  examines  the  boundary  layer 
inflow  from  the  environment  (class  II)  by  considering  the  properties  of  trajectories  that 
spend  some  time  in  the  eye  vs  those  that  are  transported  immediately  into  the  eyewall. 
Section  3c  examines  the  persistence  of  air  in  the  mid-  and  upper-tropospheric  eye  (class 
III)  by  counting  the  number  of  trajectories  that  stir  into  the  eyewall  vs  those  that  remain 
in  the  eye.  Section  3d  examines  the  ability  of  mid-  and  upper-tropospheric  environmental 
air  to  stir  into  the  eye  or  eyewall  (class  IV). 

Model  physics  were  similar  to  the  coarse-grid  runs  ex¬ 
cept  that  no  cumulus  scheme  was  used.  Further  details 
of  the  model  setup,  initialization,  and  model  physics  are 
provided  in  Braun  et  al.  (2006).  For  this  paper,  t  =  0  h 
will  be  considered  the  start  of  the  fine  mesh  simulation; 
that  is,  1800  UTC  22  August  1998. 

This  study  focuses  on  a  5-h  period  from  t  =  15-20  h, 

0900  to  1400  UTC  23  August  1998.  The  simulated  storm 
during  this  time  undergoes  an  intensifying  trend  (Fig. 

2).  The  minimum  central  surface  pressure  psfc  decreases 
steadily  from  961  mb  at  14  h  to  952  mb  at  20  h.  Maxima 
of  tangential  wind  us(c  at  z  =  40  m  during  this  time  vary 
between  50  and  55  m  s-1,  although  there  is  a  slight 


increasing  trend  in  the  6-h  period.  The  maximum  azi- 
muthally  averaged  tangential  wind  usfc  increases  slightly 
from  approximately  40  to  44  ms-1.  Comparison  with 
National  Hurricane  Center  (NHC)  best  track  for  Hur¬ 
ricane  Bonnie  at  0600,  1200,  and  1800  UTC  shows  psfc 
intensifying  (960,  958,  955  mb)  and  wind  speeds  increas¬ 
ing  (49,  51,  51  m  s-1). 

b.  Summary  of  simulated  structure 

A  complete  summary  of  the  simulated  structure  of 
the  Bonnie  simulation  can  be  found  in  Braun  et  al. 
(2006).  Several  features  of  particular  interest  to  the  cur¬ 
rent  study,  however,  are  summarized  below. 


June  2007 


CRAM  ET  AL. 


1839 


Fig.  2.  Time  series  of  minimum  central  surface  pressure  (thick  line),  maximum  tangen¬ 
tial  wind  (thin  line),  and  maximum  azimuthally  averaged  tangential  wind  (dotted  line)  at 
the  lowest  model  level  (z  =  40  m)  for  the  period  within  the  Hurricane  Bonnie  (1998) 
simulation  of  Braun  et  al.  (2006)  used  here.  The  period  covered  here  corresponds  to 
0800-1400  UTC  23  Aug  1998. 


The  system-scale  structure  of  the  storm  is  primarily 
influenced  by  northwesterly  environmental  deep  layer 
(850-200  mb)  vertical  shear.  Because  vertical  wind 
shear  has  been  identified  (DeMaria  and  Kaplan  1994) 
as  a  leading  predictor  of  storm  intensity  and  is  statisti¬ 
cally  associated  with  weaker  storms,  we  anticipate  that 
shear  will  be  an  important  factor  in  determining  simu¬ 
lated  intensity.  Simulated  radar  reflectivity  (shown  in 
Fig.  6  in  Braun  et  al.  2006)  exhibits  an  asymmetric 
storm  structure  consistent  with  the  influence  of  north¬ 
westerly  vertical  wind  shear.  The  highest  reflectivities 
appear  southeast,  east,  and  north  of  the  storm  center, 
which  is  downshear  and  left  of  the  850-200-mb  vertical 
shear  vector  calculated  about  the  storm  center.  The 
vertical  wind  shear  was  calculated  inside  a  radius  of  300 
km  from  the  storm  center  on  the  6-km  horizontal  grid 
and  has  a  magnitude  of  approximately  12  m  s-1  be¬ 
tween  z  =  1  and  z  =  12  km;  observations  from  other 
TCs  have  shown  that  at  this  magnitude,  vertical  shear 
effects  dominate  subsequent  storm  evolution  (e.g.. 
Black  et  al.  2002).  The  boundary  layer  is  characterized 
by  inflow  (Fig.  3a)  in  all  quadrants  outside  the  radius  of 
maximum  wind  (RMW),  but  is  stronger  southeast,  east, 
and  north  of  the  vortex  center.3  Generally,  the  peak 
near-surface  inflow  is  confined  below  z  =  1  km,  de¬ 
creasing  to  z  =  750  m  near  the  RMW.  In  all  regions  of 


3  The  method  for  computing  the  storm  center  is  described  in 
appendix  A. 


the  low-level  eye  except  north  of  the  center,  radial  out¬ 
flow  prevails  and  reaches  magnitudes  in  excess  of 
5  ms-1.  Radial  flow  into  and  out  of  the  core  of  the 
storm  can  be  found  at  other  levels  correlated  with  the 
expected  wave-1  pattern  associated  with  the  environ¬ 
mental  mean  flow  profile  that  contributes  strongly  to 
the  vertical  wind  shear. 

Just  above  the  surface  inflow  layer  (z  =  1.1  km),  the 
region  outside  the  RMW  exhibits  radial  outflow  to  the 
west  and  inflow  to  the  east  of  the  storm  center.  Inside 
the  RMW,  radial  outflow  exists  directly  downshear  and 
left  of  the  deep  layer  mean  vertical  shear  vector.  The 
outflow  inside  the  RMW  at  these  two  levels  (z  =  453  m 
and  1.1  km)  is  associated  with  the  convective  updrafts 
which  preferentially  form  in  the  downshear  side  of  the 
eyewall,  as  described  in  Braun  et  al.  (2006). 

In  the  middle  and  upper  troposphere  (z  =  5.6  and  9.9 
km)  radial  outflow  occurs  south,  east,  and  northeast  of 
the  storm  center.  Inside  the  RMW,  the  outflow  is  par¬ 
ticularly  strong  ( u  >  10  ms"1)  directly  south  of  the 
storm  center.  This  region  was  shown  in  Braun  et  al. 
(2006)  to  be  a  favorable  area  for  the  initiation  of  eye- 
wall  convective  hot  towers.  While  the  updrafts  presum¬ 
ably  form  at  low  levels  (below  2  km),  Braun  et  al. 
(2006)  noted  that  these  eyewall  updrafts  extended 
through  the  depth  of  the  troposphere  shortly  after  their 
initiation. 

The  azimuthal  mean  storm  structure  (Fig.  4)  gener¬ 
ally  exhibits  an  outward  sloping  eyewall  with  maximum 
mean  tangential  wind  v  of  55  m  s_1  and  mean  vertical 
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Fig.  3.  Contours  of  radial  wind  at  z  =  (a)  453  m,  (b)  1.1  km,  (c)  5.6  km,  and  (d)  9.9  km  at  t  =  15  h.  Positive  values  are  shaded  at 
levels  of  0,  2.5,  5,  and  10  m  s_1.  Negative  values  are  contoured  at  —15,  —10,  -5,  and  —2.5  m  s_1.  The  dotted  line  represents  the  radius 
of  maximum  mean  tangential  wind  at  each  level. 


velocity  vv  of  1.2  ms-1.  The  RMW  near  the  surface  at 
this  time  (f  =  15  h)  is  44  km  and  gradually  moves  in¬ 
ward  with  time  to  38  km  at  t  —  20  h  (not  shown). 
Maximum  w  increases  from  1.2  to  nearly  2.3  ms-1 
through  the  period  t  =  15-20  h. 

The  distribution  of  0e  air  (Fig.  4c,  following  Bolton 
1980)  suggests  a  reservoir  of  high  de  in  the  boundary 
layer  of  the  eye  that  is  characteristic  of  mature  hurri¬ 
canes  (W98;  Zhang  et  al.  2002).  In  the  eye,  the  maxi¬ 
mum  in  de  is  greater  than  376  K  near  the  surface  and 
reaches  a  minimum  at  z  =  4  km. 


3.  Use  of  trajectories  to  study  transport,  stirring, 
and  mixing  processes  within  the  vortex 

Each  of  the  following  sections  will  illustrate  the  in¬ 
teraction  between  various  regions  of  the  model  hurri¬ 
cane  using  distinct  sampling  strategies  for  the  initial 
locations  of  trajectories.  The  strategies  are  detailed 
separately  in  each  section.  A  simplified  view  of  a  mod¬ 
erately  sheared  hurricane  is  shown  in  Fig.  1  and  also 
describes  the  transport  processes  that  we  are  interested 
in  analyzing  in  this  study.  The  calculation  of  trajectories 
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Fig.  4.  Azimuthal  averages  of  (a)  tangential  wind,  (b)  vertical 
velocity,  and  (c)  equivalent  potential  temperature  at  t  =  15  h  in 
the  Hurricane  Bonnie  simulation  by  Braun  et  al.  (2006).  Gray 
shading  denotes  azimuthally  averaged  relative  humidity  greater 
than  or  equal  to  90%.  The  dashed  line  is  the  radius  of  maximum 
azimuthally  averaged  tangential  wind  as  a  function  of  height. 

is  performed  by  postprocessing  the  model  output  data, 
which  has  a  temporal  output  increment  of  three  min¬ 
utes.  Trajectory  positions  are  computed  every  10  s  by 
the  technique  described  in  appendix  B. 

a.  The  low-level  eye-eyewall  interaction 

1)  The  TRAJECTORY  SAMPLE 

In  this  section,  we  examine  the  interaction  between 
the  low  level  eye  and  the  eyewall  using  forward  trajec¬ 


tories  (class  I  from  Fig.  1)  seeded  in  the  moist  air  of  the 
low-level  eye  (below  z  =  1.1  km).  The  trajectories  are 
seeded  at  z  =  453  m  in  the  eye  at  t  =  900  min  =  15  h 
from  the  start  of  the  inner-grid  simulation  and  are  cal¬ 
culated  forward  for  5  h  from  this  time.  The  starting 
locations  are  on  a  rectilinear  grid  with  a  2-km  spacing  in 
both  x  and  y,  matching  the  specification  of  the  model 
grid. 

The  eye  is  distinguished  from  the  eyewall  using  the 
following  axisymmetric  criteria.  The  eyewall  is  defined 
here  as  the  region  in  the  vicinity  of  the  RMW  with 
azimuthally  averaged  vertical  velocity  W  >  0.25  ms-1 
and  axisymmetric  mean  relative  humidity  lK  >  90%. 
The  low-level  eye  (below  z  =  1  km)  is  defined  as  the 
region  inside  the  RMW  and  near  the  storm  center  with 
0e  >  370  K.  The  eye  at  z  =  1.1  km  is  defined  as  the 
region  inside  the  RMW  and  near  the  storm  center  with 
w  <  0.25  m  s-1.  The  eye  at  mid-  to  upper  levels  (z  =  5.5 
km  and  above)  is  defined  as  the  region  inside  the  RMW 
and  near  the  storm  center  with  w  <  0.25  m  s-1  and  H < 
90%.  These  definitions  using  axisymmetric  properties 
help  distinguish  the  eyewall  from  short-lived  cumuli  in 
the  eye  itself.  Although  it  would  be  preferable  to  apply 
the  same  definition  for  the  eye  at  z  =  1.1  km  as  is  used 
at  mid-  to  upper  levels,  the  relatively  higher  values  of 
at  1.1  km  (~90%-92%)  do  not  provide  a  clear  distinc¬ 
tion  between  the  mean  eye  and  eyewall  and  thus  pre¬ 
cludes  the  use  of  a  humidity  discriminant  at  this  level. 
The  shear-induced  asymmetry  in  the  storm  lowers  the 
axisymmetric  relative  humidity  below  saturation,  hence 
the  value  of  th(  =  90%  is  chosen  as  the  humidity  dis¬ 
criminant  for  defining  the  mean  eyewall  (cf.  Fig.  4). 

2)  Transport  from  the  low-level  eye  to 

THE  EYEWALL 

A  small  sample  of  the  trajectories  initiated  in  the  eye 
that  enter  the  intense  eyewall  updraft  are  shown  in  Fig. 
5.  After  first  encountering  the  eyewall,  these  seven  tra¬ 
jectories  encounter  a  downdraft  at  z  ~  1-2  km  and 
approach  the  center  of  the  storm  before  becoming  en¬ 
trained  into  the  eyewall  updraft.  At  upper  levels  (z  a  7 
km),  six  of  these  trajectories  detrain  into  the  eye.  Ver¬ 
tical  velocities  of  the  eyewall  updraft  encountered  by 
these  trajectories  range  from  1-2  m  s  1  below  z  =  3  km 
to  as  much  as  10  m  s-1  above  this  height  (not  shown). 
The  initial  6e  associated  with  the  trajectories  is  between 
372  and  376  K,  confirmation  that  these  trajectories  are 
originating  from  the  reservoir  of  relatively  high  0e  lo¬ 
cated  in  the  low-level  eye  ( de  greater  than  370  K  extend 
up  to  z  =  750  m  at  the  storm  center;  cf.  Fig.  4c).  As  the 
trajectories  rise  in  the  eyewall  updraft,  the  associated  0e 
decreases  approximately  5  K.  The  vertical  profile  of  0e 
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Fig.  5.  Seven  sample  trajectories  of  class  I  (Fig.  1)  viewed  in  (a) 
the  horizontal  x-y  plane  and  (b)  the  r-z  plane.  Trajectories  are 
seeded  less  than  15  km  from  the  storm  center  and  at  a  height  of 
453  m.  Trajectory  seed  time  is  15  h  into  the  2-km  grid  simulation 
and  trajectories  are  calculated  5  h  forward  in  time,  (c)  Associated 
values  of  Qe  as  a  function  of  height  for  the  trajectories.  Arrows 
indicate  the  direction  of  motion  of  the  trajectories.  The  5-h  time 
average  of  the  radius  of  maximum  azimuthally  averaged  tangen¬ 
tial  wind  as  a  function  of  height  is  plotted  as  a  dashed  line  in  (b). 
Seed  locations  are  marked  by  light  gray  circles. 


along  the  trajectories  becomes  relatively  constant  with 
height  above  z  =  2  km.  These  parcels  of  air,  originating 
in  the  low-level  eye,  retain  a  degree  of  warmth  relative 
to  other  eyewall  parcels,  as  is  seen  by  plotting  asym¬ 
metric  0e  along  the  trajectories  (Fig.  6a).  Asymmetries 
in  this  paper  are  defined  by  subtracting  the  axisymmet- 
ric  mean  <),,  field  from  the  complete  Qe  field.  Below  z  = 
3  km,  asymmetric  w  (Fig.  6b)  is  both  negative  and  posi¬ 
tive,  the  former  being  associated  with  the  downward 
and  inward  motion  exhibited  in  Fig.  5b.  Asymmetric  w 
is  positive  between  z  =  3  km  and  7  km,  with  values  as 
high  as  10  m  s-1.  The  asymmetric  pictures  of  0e  and  w 
suggest  that  air  stirred  into  the  eyewall  from  the  low- 
level  eye  retains  some  degree  of  positive  buoyancy  rela¬ 
tive  to  other  default  parcels  of  eyewall  air  that  did  not 
previously  reside  in  the  low-level  eye.  This  relative 
buoyancy  can  be  associated  with  the  hot  towers  identi¬ 
fied  by  Braun  (2002)  and  Smith  et  al.  (2005). 

Consider  now  the  complete  subset  of  trajectories  that 
make  their  way  out  to  the  axisymmetric  eyewall  (Fig. 
7).  There  is  an  overall  tendency  for  these  trajectories  to 
be  stirred  out  into  the  eyewall  downshear  and  left  of  the 
shear  vector,  and  stronger  tendency  is  shown  for  those 
trajectories  (dark  dots)  that  are  stirred  out  near  the 
base  of  the  eyewall.  Of  the  551  trajectories  that  are 
seeded  inside  the  eye  at  this  level,  313  trajectories,  or 
56.8%,  are  mixed  out  to  the  eyewall  within  5  h.  Al¬ 
though  the  escape  points  shown  are  located  all  around 
the  RMW  azimuth  ring,  the  majority  of  trajectories 
(Fig.  7b)  encounter  the  eyewall  at  locations  southeast, 
east,  and  northeast  of  the  storm  center.  These  findings 
are  consistent  with  the  results  in  Braun  et  al.  (2006), 
which  show  that  the  main  eyewall  updrafts  occur  down- 
shear  and  left  of  the  mean  layer  vertical  shear  vector. 

3)  Evidence  of  superintensity  mechanism 

PM03  defined  superintensity  as  a  state  where  hurri¬ 
cane  intensity  exceeds  the  predicted  maximum  inten¬ 
sity  of  Emanuel’s  (1995)  maximum  potential  intensity 
(E-PI)  theory.  The  mechanism  identified  by  PM03  in 
their  axisymmetric  modeling  was  transport  of  air  from 
the  low-level  eye  to  the  eyewall.  The  low-level  eye  air 
was  enhanced  in  0e  relative  to  the  eyewall  air;  thus,  its 
introduction  to  the  base  of  the  eyewall  represented  an 
additional  source  of  heat  for  the  hurricane  heat  engine 
cycle.  In  E-PI  theory,  the  only  source  of  heat  for  driving 
the  hurricane  circulation  is  the  ocean  beneath  the  eye- 
wall  through  a  wind-induced  vertical  heat  exchange. 

The  representative  sample  of  trajectories  shown  in 
Fig.  5  provides  evidence  of  the  superintensity  mecha¬ 
nism  operating  in  the  MM5  Bonnie  simulation.  Once 
these  sample  trajectories  are  stirred  out  to  the  eyewall, 
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Fig.  6.  (a)  Asymmetric  8e  and  (b)  asymmetric  vt>  as  a  function  of  height  along  the  trajectories  of  class 

I  (Fig.  1)  shown  in  Fig.  5. 


they  are  lifted  by  the  eyewall  updrafts.  The  form  of  the 
interaction  between  the  relatively  high-0c  trajectory 
particles  and  the  eyewall  air  is  an  incomplete  reduction 
of  trajectory  0e  to  the  mean  eyewall  value,  retaining  a 
moderate  localized  enhancement  of  the  eyewall  updraft 
associated  with  the  trajectory  (Fig.  6).  Since  0e  is  ap¬ 
proximately  conserved,4  a  reduction  of  0e  following  the 
trajectory  upon  introduction  to  the  eyewall  must  be 
associated  with  an  increase  in  0e  in  the  surrounding  air, 
thus  warming  the  axisymmetric  mean  0e  by  mixing. 

The  stirring  of  high-0e  air  from  the  low-level  eye  into 
the  eyewall  is  corroborated  by  instantaneous  r—z  cross 
sections  of  0e  along  radial  legs  in  the  southeast  quadrant 
of  the  storm  (not  shown).  As  shown  in  Braun  et  al. 
(2006),  the  eyewall  updrafts  originate  in  this  quadrant 
of  the  storm,  which  corresponds  to  the  downtilt-right- 
hand  side  of  the  eyewall  in  relation  to  the  environmen¬ 
tal  vertical  shear  vector.  These  cross  sections  of  0e  ex¬ 
hibit  a  tongue  of  high-0e  air  in  the  low-level  eye  extend¬ 
ing  into  the  base  of  the  eyewall  updrafts  [see  an 


4  The  MM5  model  is  not  designed  to  conserve  8e  of  Bolton 
(1980).  Both  8e  and  MM5,  however,  are  designed  in  reference  to 
reality,  so  approximate  conservation  can  be  anticipated  (and  is 
roughly  observed  in  the  post  analyses)  in  the  absence  of  ice  pro¬ 
cesses,  surface  interaction,  or  strong  subgrid  scale  diffusive  mix¬ 
ing. 


example  of  this  structure  in  Fig.  12  of  Braun  (2002)]. 
This  suggests  that  this  warmer  0e  air  is  being  introduced 
to  the  eyewall  updrafts  as  they  form;  it  was  not  deter¬ 
mined  whether  the  introduction  of  the  high  0,,  air  into 
the  eyewall  was  the  cause  of  the  updraft  initiation. 

For  the  case  of  the  Bonnie  simulation,  the  sea  surface 
temperature  is  303  K,  the  latitude  is  24.2°N,  and  the 
ambient  surface  pressure  is  1010  mb.  Since  there  is  no 
dissipative  heating  in  the  numerical  simulation, 
Emanuel  (1995)  represents  the  appropriate  derivation 
of  E-PI  for  comparison.  A  reasonable  outflow  tempera¬ 
ture  following  the  technique  of  PM03  of  the  axisym¬ 
metric  circulation  is  210  K.  To  supply  the  values  of 
surface  relative  humidity,  we  recall  the  closure  assump¬ 
tion  invoked  by  Emanuel  (1995)  that  the  relative  hu¬ 
midity  is  hypothesized  to  have  the  same  value  at  the 
eyewall  as  in  the  environment.  We  choose,  therefore, 
9{  —  80%  as  an  expedient  representative  of  the  envi¬ 
ronment,  in  the  spirit  of  the  a  priori  approach  of  E-PI. 
Figure  4  shows,  however,  that  relative  humidity  can  be 
greater  than  90%  in  the  subcloud  inflow  layer.  We  thus 
also  provide  here  a  calculation  using  “H  =  90%  as  rep¬ 
resentative  of  the  eyewall,  but  we  can  anticipate  a 
weaker  magnitude  of  E-PI  for  a  higher  H  (Camp  and 
Montgomery  2001).  Diagnosis  of  CkICD  for  the  surface 
flux  parameterization  for  this  model  run  exhibits  a  spa¬ 
tial  variability  of  the  ratio,  ranging  from  0.65  for  the 
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Fig.  7.  (a)  Locations  where  trajectories  seeded  inside  the  eye  ( de  >  370  K)  at  a  height  of 
z  =  453  m  (class  I  of  Fig.  1)  are  stirred  out  to  the  axisymmetric  eyewall.  The  arrow  in  the  upper 
right  comer  indicates  the  direction  of  the  mean  layer  (850-200  mb)  vertical  shear  vector. 
Eyewall  encounters  occurring  at  or  below  z  =  1.5  km  are  denoted  by  dark  circles,  encounters 
above  z  =  1.5  km  are  light  gray,  (b)  Histogram  of  the  azimuths  of  the  points  shown  in  (a).  Of 
the  551  total  trajectories  seeded  inside  the  eye  at  this  level,  313  trajectories  (57%)  are  stirred 
out  and  meet  the  eyewall  criteria.  Eye  seed  points  extend  out  to  a  radius  of  24  km. 


average  over  the  complete  inner  domain  to  0.35  at  the 
radius  of  maximum  winds.  With  these  values  and  using 
tH  =  80%,  E-PI  predicts  values  of  maximum  v  at  the 
top  of  the  boundary  layer  ranging  from  52.2  m  s-1 
for  CkICD  =  0.65  to  38.2  m  s-1  for  Ck/CD  =  0.35.  For 
=  90%,  the  values  are  36.5  m  s_1  for  CkICD  =  0.65 
and  26.9  m  s-1  for  CkICD  =  0.35.  The  strongest  axisym¬ 
metric  tangential  winds  of  55  m  s-1  found  at  ( r  =  48, 
z  =  0.8  km)  in  the  Bonnie  simulation  are  near  to,  or 
much  greater  than,  E-PI.  The  simulation  presented 
here  is  thus  plausibly  superintense.  Considering  that 
the  simulated  Bonnie  presented  here  was  under  the 
influence  of  moderate  vertical  wind  shear  [which  is 
associated  with  weaker  intensity  in  the  literature  (De- 
Maria  1996;  Frank  and  Ritchie  2001;  Wu  and  Braun 
2004)],  it  is  noteworthy  that  the  storm  is  near  its  E-PI 
value  (or  much  stronger,  depending  on  selection  of  pa¬ 
rameters).5 


5  In  a  strict  sense,  E-PI  theory  is  valid  for  tropical  cyclones  in  a 
steady-state  configuration.  Although  the  maximum  tangential 
wind  speed  exhibits  a  near-steady  tendency  (Fig.  2),  the  minimum 
central  pressure  decreases  nearly  linearly  with  time,  thus  suggest¬ 
ing  the  simulated  storm  is  not  in  a  steady-state  configuration  dur¬ 
ing  our  analysis  period.  Given  the  thermodynamic  conditions,  one 
can  nevertheless  compare  the  modeled  storm  intensity  to  the  com¬ 
parable  equilibrium  intensity  anticipated  by  theory.  While  we  ac¬ 
cept  that  transients  could  increase  intensity  stronger  than  E-PI  for 
some  short  period  of  time,  we  suggest  that  only  the  eye-eyewall 
exchanges  of  warmth  and  moisture  can  sustain  the  storm  at  or 
stronger  than  E-PI  for  such  a  multihour  period  in  the  face  of 
moderate  vertical  wind  shear. 


b.  Eye  inflow  trajectories  and  replenishment  of  the 
6e  reservoir 

1)  The  TRAJECTORY  SAMPLE 

If  the  eyewall  updrafts  are  consuming  some  of  the 
high-0,,  air  within  the  low-level  eye,  then  there  must  be 
some  mechanism  for  replenishing  that  high-0,,  reservoir. 
To  examine  this  replenishment  process,  we  investigate 
forward  trajectories  (class  II  of  Fig.  1)  seeded  in  the  in¬ 
flow  layer  converging  on  the  storm  center.  The  question 
to  be  addressed  at  this  point  is  what  percentage  of  the 
trajectories  is  transported  underneath  the  eyewall  up¬ 
drafts  and  into  the  hurricane  eye,  as  opposed  to  following 
the  standard  ascent  path  through  the  eyewall  updrafts? 

Seed  points  are  located  horizontally  on  a  radius  ver¬ 
sus  azimuth  grid.  Radial  seed  spacing  is  5  km,  between 
50  <  r  <  180  km.  Twenty-five  azimuthal  seed  points 
exist  at  each  seed  radius.  The  vertical  seed  levels  used 
are  z  =  40, 121,  and  244  m.  This  seeding  strategy  results 
in  625  trajectories  calculated  at  each  vertical  seed  level. 
The  trajectories  are  calculated  forward  3  h  from  a  seed 
time  of  t  =  15  h  into  the  2-km  simulation. 

To  ensure  that  trajectories  seeded  at  or  near  the  low¬ 
est  model  level  ( z  =  40  m)  do  not  cross  the  lower  grid 
boundary,  it  was  necessary  for  us  to  extrapolate  the 
velocity  data  to  a  near-surface  height;  we  chose  z  =  10 
m.  Here  we  follow  the  logarithmic  boundary  layer  as 
observed  by  Powell  et  al.  (2003)6  and  assign  the  mag- 


6  We  assume  that  the  parameterization  of  the  surface  layer  in 
the  model  follows  similar  assumptions. 
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Fig.  8.  Histogram  of  trajectories  of  class  II  (Fig.  1)  that  are  transported  into  the  eye,  shown  as  a 
function  of  trajectory  seed  radius.  Trajectory  seed  points  for  this  analysis  are  every  5  km  in  the  radial 
direction  and  seed  levels  of  z  =  40,  121,  and  244  m.  Bin  size  is  5  km  and  sample  size  is  75  trajectories 
in  each  bin  (25  at  each  seed  level).  Percentages  are  shown  for  (a)  seed  level  z  =  244  m,  (b)  seed  level 
z  =  121  m,  and  (c)  seed  level  z  =  40  m.  A  trajectory  is  considered  to  be  transported  into  the  eye  if  its 
trajectory  8e  >  370  K  at  z  <  2  km. 


nitudes  of  the  u-  and  i>velocity  components  to  be  85% 
of  their  values  at  z  =  40  m.  We  assign  the  vertical 
velocity  w  to  vanish  at  z  =  10  m  in  order  to  ensure  that 
no  trajectories  descend  below  10  m. 

2)  Low-level  eye  mass  replacement 

For  a  trajectory  of  class  II  to  be  transported  into  the 
eye,  we  required  the  0e  value  associated  with  a  particu¬ 
lar  trajectory  to  exceed  370  K  below  a  vertical  height  of 
z  —  2  km.  This  threshold  value  is  chosen  based  on  the 
values  of  de  in  the  high-0,,  reservoir  of  the  low  level  eye 
(Fig.  4c).  By  this  criterion,  a  significant  fraction  of  class 
II  trajectories  (i.e.,  the  boundary  layer  inflow)  spends  at 
least  some  time  in  the  eye  (Fig.  8).  Considering  trajec¬ 
tories  originating  just  outside  the  eyewall  ( r  <  100  km), 
close  to  half  of  low  level  (z  =  40  m  seeding;  Fig.  8c)  air 


is  found  to  slip  under  the  eyewall  and  enter  the  eye, 
with  probability  decreasing  with  height  to  about  20%  at 
z  =  121  m.  At  further  radii  (r  >  125  km),  turbulent 
redistribution  between  levels  in  the  boundary  layer 
makes  the  distinction  between  the  percentages  at  these 
levels  less  definitive.  Additionally,  at  larger  radii,  there 
is  an  increased  chance  that  trajectories  will  be  intercepted 
by  the  outer  rain  bands  or  that  the  3-h  integration  is 
insufficient  to  permit  the  inflow  to  reach  the  eyewall. 

We  now  seek  trajectories  that  can  serve  as  prototypes 
for  the  two  primary  consequences  for  class  II  trajecto¬ 
ries  that  reach  the  eyewall:  either  to  enter  the  eye  for 
some  period  of  time  or  to  be  immediately  drawn  into 
the  eyewall  updraft  without  encountering  the  eye.  The 
approach  for  constructing  these  prototypes  will  elimi¬ 
nate  a  large  number  of  borderline  cases,  since  these 
trajectories  typify  neither  consequence.  For  the  eye  in- 
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Fig.  9.  Radius  vs  height  plots  for  a  select  sample  of  class  II  (Fig.  1)  trajectories  that  can  be  identified  as  (a)  eye 
reservoir  trajectories  or  (b)  standard  eyewall  ascent  trajectories,  and  (c),(d)  0e  vs  height  for  the  same  trajectories. 
Sample  sizes  are  41  for  eye  reservoir  and  78  for  standard  eyewall  ascent. 


flow  prototypes,  consider  those  trajectories  of  class  II 
defined  above  which  exceed  de  =  372  K  at  some  point 
in  the  calculation.  For  the  eyewall  ascent  prototypes, 
consider  those  trajectories  of  class  II  with  maximum  6e 
between  369  K  <  0e  <  370  K  (Fig.  9).  The  amount  of 
time  eye  inflow  prototype  trajectories  spend  in  the  eye 
can  vary  greatly,  where  this  time  is  defined  as  the  pe¬ 
riod  spent  with  de  >  370  K  just  prior  to  entrainment 
into  the  eyewall  updraft.  The  amount  of  time  required 
for  a  trajectory  to  acquire  a  substantial  increase  in  de 
(3  K)  can  be  as  short  as  15  min  (Fig.  10),  although  40-60 
min  is  typical.  The  behavior  described  here  is  sufficient 
to  explain  how  the  eye  entropy  reservoir  is  maintained; 
the  mass  to  replace  air  entrained  into  the  eyewall  can 
be  adequately  accounted  for  while  it  is  demonstrated 
that  the  6e  of  these  parcels  can  be  quickly  enhanced 
through  exchange  with  the  underlying  ocean. 

3)  Superintensity  mechanism  maintained 

The  eye  inflow  trajectories  (Figs.  9a, c)  show  clear 
evidence  of  the  superintensity  mechanism  occurring. 


All  trajectories  converge  toward  the  center  near  the 
surface  ( z  <  200  m  at  r  >  40  km  and  z  —  500  m  at 
r  <  40  km)  and  approach  to  within  15-25  km  of  the 
storm  center  as  they  are  stirred  into  the  eye.  When 
they  are  stirred  out  to  the  eyewall  and  lifted  by  the 
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Fig.  10.  Maximum  t)e  vs  eye  residence  time  for  the  eye  reservoir 
trajectories  shown  in  Figs.  9a, c.  These  trajectories  are  of  class  II 
from  Fig.  1. 
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Fig.  11.  Difference  in  de  between  z  =  1.5  km  and  500  m  [A0C  =  de(z  =  1.5  km)  -  6e(z  =  0.5  km)]  for 
the  eye  reservoir  (triangles)  and  standard  eyewall  ascent  (squares)  trajectories  shown  in  Figs.  9a,c.  The 
numbering  of  trajectories  along  the  x  axis  is  arbitrary.  These  trajectories  are  of  class  II  from  Fig.  1. 


eyewall  updrafts  (Fig.  9c),  the  interaction  between 
the  high-0,,  trajectory  particles  and  the  eyewall  air  is 
apparent.  The  6e  value  of  the  trajectories  is  not  main¬ 
tained  at  the  values  found  in  the  low-level  eye,  but 
decrease  partially  to  the  value  of  the  surrounding  eye- 
wall  air  as  a  result  of  diffusive  mixing.  Conversely,  the 
eyewall  ascent  trajectories  (Figs.  9b,d)  show  a  slight 
warming  initially  in  0e  as  they  ascend  in  the  eyewall 
updraft. 

After  exiting  the  eye,  most  of  the  eye  inflow  trajec¬ 
tories  exhibit  6e  decreases  of  as  much  as  4  K  (Fig.  11), 
as  a  result  of  mixing  with  relatively  cooler  eyewall  air. 
The  eyewall  ascent  trajectories,  in  turn,  exhibit  in¬ 
creases  in  6e — up  to  3  K — resulting  from  their  interac¬ 
tion  with  the  higher- 0e  air  coming  from  the  low-level 
eye. 

As  the  moist  entropy  of  the  low-level  eye  is  stirred 
into  to  the  eyewall,  the  reservoir  of  high-0,,  air  in  the 
low-level  eye  must  be  replenished.  This  is  achieved  pri¬ 
marily  via  the  convergence  of  air  toward  the  storm  cen¬ 
ter  in  the  low-level  inflow  layer  as  shown  above  (Fig. 
9a). 


c.  The  upper-level  eye-eyewall  interaction 

1)  The  trajectory  sample 

To  study  the  degree  of  isolation  of  the  eye  from  the 
eyewall  in  the  mid-  to  upper  troposphere,  we  examine 
in  this  section  forward  trajectories  (class  III  of  Fig.  1) 
seeded  within  the  eye  above  the  inversion.  The  proce¬ 
dure  follows  section  3a  except  that  seeding  at  heights  of 
1.1,  5.6,  and  9.9  km  will  also  be  considered.  The  latter 
two  are  taken  to  represent  air  parcels  from  the  mid-  to 


upper-level  eye,  while  1.1  km  is  more  representative  of 
the  lower-level  eye. 

2)  The  containment  vessel  hypothesis 

REVIEWED 

For  trajectories  originating  inside  the  eye  in  the 
middle  and  upper  troposphere,  we  are  interested  in 
whether  the  eye  behaves  as  a  containment  vessel.  Pro¬ 
posed  by  W98,  the  containment  vessel  hypothesis  posits 
that  the  air  contained  inside  the  eye  above  the  eye  in¬ 
version  in  intense  tropical  cyclones  has  remained  inside 
the  eye  since  the  eye  formed.  Using  flight-level  and 
dropsonde  data,  W98  identified  eye  soundings  typical 
of  intensifying  tropical  cyclones.  The  main  features  of 
these  soundings  included  an  eye  inversion  between  900 
and  850  mb  and  dewpoint  depressions  of  approximately 
10  K  above  the  inversion.  The  rise  and/or  fall  of  the 
inversion  results  from  a  balance  between  low-level 
moist  inflow,  which  is  brought  into  the  eye  to  compen¬ 
sate  for  the  moist  low-level  eye  air  lost  to  the  eyewall 
updrafts  (as  shown  previously  in  section  3b),  and  warm, 
dry  subsidence  above  the  inversion  induced  by  the  loss 
of  low-level  mass  to  the  eyewall  updrafts.  W98  de¬ 
scribes  the  convective  eyewall  updrafts  as  heat  pumps, 
which  do  work  on  the  eye  by  drawing  moist  air  out  of 
the  low-level  eye  and  hence  force  thermally  indirect 
descent  above  the  eye  inversion. 

The  eye  soundings  typical  of  weakening  hurricanes 
(W98)  showed  eyewall  moisture  being  mixed  into  the 
eye  at  mid-  to  upper  levels.  In  contrast  to  intensifying 
storms,  the  weakening  case  is  characterized  by  a  more 
gradual  Y-shaped  inversion.  The  base  of  the  inversion 
in  the  weakening  cases  is  typically  found  (W98)  be¬ 
tween  850  and  700  mb. 
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Fig.  12.  Skew  T-log  p  diagram  of  the  sounding  located  inside 
the  eye  at  the  storm  center.  The  time  displayed  is  t  =  17.5  h,  or 
halfway  through  the  5-h  forward  trajectory  calculation.  Tempera¬ 
ture  is  the  solid  line  and  dewpoint  is  the  dashed  line. 


3)  Evidence  of  upper-level  eye  mass 

RECYCLING 

The  MM5  Bonnie  simulation  described  in  Braun  et 
al.  (2006)  and  herein  appears  to  be  more  typical  of  the 
soundings  in  W98  in  which  eyewall  moisture  is  being 
stirred  into  the  eye  (Fig.  12).  The  center  eye  sounding 
follows  a  saturated  adiabat  between  900  and  650  mb, 
with  the  exception  of  a  shallow  quasi-isothermal  layer 
between  850  and  800  mb.  The  air  between  850  mb  and 
the  ocean  is  moister,  with  a  fairly  constant  dewpoint 
depression  of  approximately  0.8  K.  Above  the  isother¬ 
mal  layer,  dewpoint  depressions  range  between  2  and 
8  K. 

Trajectories  seeded  in  the  eye  at  the  z  =  5.6-  and 
9.8-km  levels  suggest  that  more  than  half  (59%  at  5.6 
km,  76%  at  9.8  km)  of  the  eye  air  is  removed  from  the 
axisymmetric  eye  to  the  axisymmetric  eyewall  (using 
the  same  criteria  as  above)  over  the  5-h  period  of  study. 
At  z  =  5.6  km,  the  exit  locations  for  trajectories  (Fig. 
13)  tend  to  be  toward  the  southeast  sector  of  the  eye, 
consistent  with  the  radial  outflow  shown  at  these  azi¬ 
muths  (Fig.  3c).  At  z  =  9.8  km,  the  pattern  is  the  same 
but  is  shifted  toward  southern  azimuths  (Figs.  14,  3d). 
Braun  et  al.  (2006)  has  associated  this  radial  wind  sig¬ 
nature  with  the  outflow  from  strong  updraft  towers 
which  form  on  the  southern  side  of  the  eyewall. 

Figure  15  shows  the  percentage  of  eye  trajectories 
seeded  at  low-  to  midlevels  (z  =  453  m,  1.1  km,  and  5.6 
km)  that  are  stirred  into  the  eyewall.  In  general,  the 
percentage  increases  with  increasing  seed  radius,  and 


x  (km)  Direction  from  center 

Fig.  13.  (a)  Locations  where  trajectories  seeded  inside  the  eye  (w  <  0.25  m  s_1  and  7f< 
90%)  at  a  height  of  z  =  5.6  km  are  stirred  into  the  axisymmetric  eyewall.  The  arrow  in  the 
upper  right  corner  indicates  the  direction  of  the  mean  layer  (850-200  mb)  vertical  shear 
vector,  (b)  Histogram  of  the  azimuths  of  the  points  shown  in  (a).  Of  the  1497  total 
trajectories  seeded  inside  the  eye  at  this  level,  878  trajectories  (59%)  are  stirred  out  and 
meet  the  eyewall  criteria.  These  trajectories  are  of  class  III  from  Fig.  1.  Eye  seed  points 
extend  out  to  a  radius  of  40  km. 
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Fig.  14.  As  in  Fig.  13  but  at  a  height  of  z  =  9.8  km.  Of  the  1601  total  trajectories  seeded 
inside  the  eye  at  this  level,  1210  trajectories  (76%)  are  stirred  out  and  meet  the  eyewall 
criteria.  Eye  seed  points  extend  out  to  a  radius  of  42  km. 


the  percentage  of  trajectories  seeded  at  lower  levels 
(z  =  453  m  and  1.1  km)  that  are  being  stirred  into  the 
eyewall  is  approximately  10%-20%  higher  than  what  is 
shown  at  the  z  =  5.6  km  seed  level.  A  large  fraction 
(45%  or  more)  of  air  from  lower  levels  (z  =  453  m  and 
1.1  km)  of  the  eye  and  within  r  =  15  km  of  the  center 
encounters  the  eye  within  5  h,  evidence  for  a  large  de¬ 
gree  of  mass  replacement  of  the  eye. 

Although  we  have  not  comprehensively  documented 
the  mid-  to  upper-level  eye  replacement  (which  would 
require  a  longer  time  interval),  the  process  is  illustrated 
in  part  by  the  sample  of  seven  class  I  trajectories  (Fig. 
5),  where  air  rising  through  the  eyewall  detrains  into 
the  eye  above  z  =  6  km. 


d.  Ventilation  of  the  TC  by  the  midlevel 

environment 

1)  The  TRAJECTORY  SAMPLE 

To  examine  the  degree  to  which  dry,  midlevel  envi¬ 
ronmental  air  intrudes  upon  the  eyewall  (ventilation) 
for  this  section,  we  investigate  backward  trajectories 
seeded  across  the  eyewall  and  determine  the  sources  of 
that  air.  All  four  classes  of  trajectories  found  in  Fig.  1 
can  in  principle  be  found  in  such  a  sample,  and  the 
objective  will  be  to  distinguish  class  IV  trajectories 
from  the  sample  and  examine  their  thermodynamic 
characteristics. 

The  seed  locations  are  at  z  =  5  km  with  equal  radial/ 


Trajectory  seed  radius  (km) 


Fig.  15.  Percentage  of  eye  trajectories  of  class  III  seeded  at  each  radius  and  at  z  =  453 
m  (solid  line),  1.1  km  (dotted  line),  and  5.6  km  (dashed  line),  which  satisfy  the  criteria  for 
being  stirred  into  the  axisymmetric  mean  eyewall  (w  >  0.25  m  s_1  and  7/  >  90%).  Per¬ 
centages  are  plotted  as  a  function  of  trajectory  seed  radius.  Sample  sizes  are  551  (453  m), 
903  (1.1  km),  and  1497  (5.6  km)  trajectories.  Radial  locations  are  binned  to  the  nearest 
1  km. 
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Fig.  16.  (a)  Locations  where  trajectories  encounter  the  axisymmetric  eyewall  (w  >  0.25  m  s_1)  plotted  in  the 
(r,  z)  plane.  Dark  circles  are  trajectories  determined  to  be  class  I  or  II  (boundary  layer  inflow  or  low-level  eye 
interaction;  see  Fig.  1);  triangles  are  class  III  (mid-  to  upper-level  eye  interaction);  and  light  circles  are  class  IV 
(mid-  to  upper-level  environmental  interaction),  (b)  The  same  is  shown  in  plan  view.  For  each  class  IV  trajectory, 
the  eyewall  entry  point  is  again  plotted  with  light  circles  and  the  point  of  farthest  radius  is  shown  with  X’s  (c)  in  an 
(r,  z)  plane  and  (d)  as  a  plan-view  plot,  with  trajectories  overlaid.  The  dashed  line  is  the  radius  of  maximum  updraft. 


azimuthal  spacing  between  36  <  r  <  50  km.  There  are 
50  seed  points  at  each  radius,  and  the  radial  spacing  is 
2  km,  for  a  total  of  400  trajectory  seeds.  The  seed  time 
used  is  t  =  20  h  and  the  backward  calculation  is  carried 
out  for  5  h.  Throughout  the  seeding  annulus  w  >  0.25 
m  s-1. 

The  process  for  distinguishing  different  classes  of  tra¬ 
jectories  is  as  follows.  For  each  back  trajectory,  w  > 
0.25  m  s-1  at  the  seed  (final)  time,  by  construction. 
Each  trajectory  may  be  followed  backward  in  time  from 
the  seed  location  to  determine  the  first  location  where 
w  is  no  longer  greater  than  >0.25  m  s-1,  which  we  de¬ 
fine  as  the  eyewall  entry  location.  It  is  possible  that  one 
or  more  trajectories  might  remain  in  the  axisymmetric 


eyewall  throughout  the  trajectory  calculation,  but  it  was 
found  that  all  trajectories  had  origins  outside  the  eye- 
wall  through  the  5-h  calculation  period.  Trajectories 
with  eyewall  entry  points  below  z  =  1.5  km  are  consid¬ 
ered  to  be  in  either  class  I  or  II,  and  no  effort  is  made 
in  the  rest  of  this  section  to  distinguish  class  I  from  class 
II.  Trajectories  with  eyewall  entry  above  z  =  1.5  km  but 
at  a  radius  interior  to  the  radius  of  maximum  updraft 
are  considered  class  III.  For  the  rest  of  the  trajectories 
(eyewall  entry  above  z  =  1.5  km  and  outside  the  radius 
of  maximum  updraft),  two  modes  were  identified  (Fig. 
16a).  The  first  mode  of  trajectories  were  involved  in 
updrafts  just  exterior  to  the  axisymmetric  eyewall  (be¬ 
cause  of  a  noncircular  shape  of  the  eyewall)  that 
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emerged  from  the  boundary  layer;  these  were  classified 
with  the  class  I/II  trajectories  by  their  starting  points 
being  below  z  =  1.5  km.  The  remainder  of  trajectories 
were  then  classified  as  class  IV;  thus  these  are  trajec¬ 
tories  with  eyewall  entry  points  above  z  =  1.5  km,  radii 
exterior  to  the  radius  of  maximum  updraft,  and  starting 
locations  above  z  =  1.5  km. 

2)  Environmental  entrainment 

Of  the  400  trajectories,  213  are  of  class  I/II,  115  are  of 
class  III,  and  72  are  of  class  IV.  Class  I/II  trajectories 
are  somewhat  evenly  distributed  in  azimuth  (dark  dots; 
Fig.  16b).  The  distribution  of  class  III  (triangles)  and 
class  IV  (light  dots;  Fig.  16b)  on  the  other  hand  show 
distinct  azimuthal  biases,  apparently  related  to  the  12 
m  s"  1  northwesterly  shear  across  the  system;  class  III 
entry  points  appear  on  the  downshear  and  class  IV  en¬ 
try  points  on  the  upshear  azimuths  of  the  storm.  The 
locations  where  the  class  IV  trajectories  have  their 
maximum  radius  can  be  used  (X  symbols  in  Figs.  16c, d) 
to  study  the  means  of  interaction  from  the  environ¬ 
ment.  These  maximum  radius  points  are  generally  po¬ 
sition  to  the  left  of  the  shear  vector.  Many  of  these 
trajectories  quickly  make  their  way  to  the  eyewall. 
within  a  quarter  rotation  about  the  center;  while  others, 
whose  maximum  radii  occur  up  to  five  full  hours  prior 
to  the  seed  time  (the  maximum  possible  by  this  calcu¬ 
lation),  orbit  more  than  2  times  before  encountering 
the  eyewall.  Several  of  those  trajectories  that  make  the 
quick  approach  to  the  eyewall  appear  to  descend  in 
approach  along  a  slanted  path,  while  the  rest  appear  to 
maintain  a  constant  height  during  approach  to  the  eye- 
wall. 

3)  Ventilation  illustrated 

Ventilation  is  thought  to  weaken  the  hurricane, 
largely  because  of  the  transport  of  properties  from  the 
environment  that  are  detrimental  to  maintenance  of  the 
hurricane.  Our  approach  here  is  to  investigate  thermo¬ 
dynamic  alterations  of  the  eyewall  by  environmental 
air.  The  midtroposphere  in  the  tropical  environment 
typically  shows  a  midlevel  (),,  minimum  (Holton  1992); 
in  terms  of  0e,  the  hurricane  eye/eyewall  are  the  warm¬ 
est  areas  outside  of  the  stratosphere.  Any  lateral  mixing 
from  the  midtropospheric  environment  into  the  eyewall 
will  reduce  eyewall  9e  on  average.  Such  a  reduction  in  6e 
has  been  hypothesized  to  reduce  the  effectiveness  of 
the  hurricane  heat  engine  (Riehl  and  Malkus  1961; 
Emanuel  et  al.  2004)  and/or  weaken  the  warm-core 
temperature  anomaly  aloft  that  supports  the  low  sur¬ 
face  pressure  signature  of  the  hurricane  (Gray  1968; 
Frank  and  Ritchie  2001;  Knaff  et  al.  2004). 


Table  1.  Midtropospheric  environmental  trajectory  (class  IV) 
statistics.  The  change  in  8e  along  the  trajectories  is  A8e  =  8e(t  = 
farthest  radius)  —  8e(t  =  seed  point).  Ensemble  averages  are  de¬ 
noted  by  (•),  cr  is  standard  deviation,  prime  symbols  denote  de¬ 
partures  from  the  axisymmetric  mean,  and  AA  represents  the 
change  in  azimuth  along  the  trajectory  paths. 


A0e  <  0 

A0„  >  0 

All 

N 

32 

40 

72 

<A0„> 

-4.87  K 

2.30  K 

-0.89  K 

4.39  K 

1.38  K 

4.72  K 

( d'e )  (seed  point) 

-1.30  K 

-5.20  K 

-3.47  K 

cr(0() 

3.79  K 

2.31  K 

3.60  K 

(A0e)  +  < 8'e > 

-6.17  K 

-2.90  K 

-4.36  K 

(AA) 

8.44  rad 

7.76  rad 

8.06  rad 

cr  (AA) 

3.29  rad 

2.58  rad 

2.92  rad 

The  simplest  measure  of  ventilation  of  eyewall  air 
with  environmental  air  is  to  note  that  72  of  400  trajec¬ 
tories,  representing  18%  of  the  mass  (and  number,  co¬ 
incidentally)  rising  through  z  =  5  km  in  the  eyewall, 
originated  in  the  environment  (light  colored  dots;  Fig. 
16).  The  large  number  of  interacting  trajectories  would 
appear  to  be  a  consequence  of  analyzing  a  sheared  hur¬ 
ricane  simulation.  While  a  direct  comparison  of  mixing 
processes  in  a  three-dimensional  simulation  and  axi- 
symmetry  must  be  done  cautiously,  we  performed  a 
comparative  axisymmetric  simulation  using  the  4x  run 
from  PM03  (which  is  necessarily  without  vertical  wind 
shear;  not  shown)  and  found  zero  class  IV  trajectories. 
We  plan,  in  future  work,  to  extend  the  present  analysis 
to  idealized,  three-dimensional  hurricane  simulations 
without  vertical  shear. 

An  analysis  of  entrant  class  IV  trajectories  is  sum¬ 
marized  in  Table  1.  Here,  the  change  in  9e  along  the 
trajectory  is  measured  using 

A0e  =  9e(t  =  farthest  radius)  —  Qe(t  =  seed  point) 

(1) 

and  the  dataset  is  partitioned  according  to  whether  A0e 
is  positive  or  negative.  When  A 9e  <  0,  the  trajectory 
warms  on  approach  to  the  eyewall,  which  must  come  at 
the  cost  of  a  reduction  of  9e  from  the  eyewall  and  the 
surrounding  region,  either  through  warming  or  moist¬ 
ening  of  the  newly  arriving  trajectory.  On  average,  class 
IV  trajectories  show  negative  values  ((A0e)  =  —0.89  K; 
where  {■)  is  used  for  an  ensemble  average),  indicating  a 
net  warming  of  such  parcels  upon  approach  to  and  en¬ 
trance  into  the  eyewall  by  some  diffusive  or  irreversible 
mixing.  Despite  this  warming,  these  class  IV  trajecto¬ 
ries  remain  cool  relative  to  the  axisymmetric  mean  (0'  < 
0).  Also  shown  in  Table  1  are  standard  deviations  cr  of 
A  0,,  and  0',.  Eyewall  9e  is  reduced  by  the  introduction  of 
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class  IV  trajectories,  either  directly  by  stirring  cold  air 
into  the  vicinity  of  the  warmer  eyewall  air  that  emerges 
from  the  boundary  layer  inflow  and  the  eye,  or  indi¬ 
rectly  by  diffusive  mixing.  Since  class  IV  represents 
about  18%  of  the  mass  at  z  =  5  km,  the  axisymmetric 
mean  eyewall  is  cooled 

0.18 

A0e(eyewall)  «=  1  Q  _  Q  lg  (-0.89  -  3.47)  K 

«  -  0.96  K.  (2) 

The  large  standard  deviations  and  the  unbalanced 
sample  size  of  positive  and  negative  A0es  are  indicative 
of  a  positively  skewed  distribution,  where,  while  all 
class  IV  parcels  stir  in  mass,  a  small  number  of  events 
may  perform  the  largest  modification  of  the  eyewall 
mean  9e. 

The  technique  of  this  section  was  applied  using  a 
second  seed  time  (not  shown),  t  —  19  h  instead  of  20  h 
used  above.  At  this  time,  the  axisymmetric  eyewall  was 
narrower,  so  the  seed  annulus  was  38  <  r  <  48  km.  This 
provides  300  trajectories  for  analysis,  81  of  which 
(27%)  were  identified  as  class  IV,  thus  there  were  more 
ventilating  air  parcels  at  this  time.  On  the  whole,  this 
set  of  trajectories  differs  from  the  set  illustrated  above 
by  a  net  decrease  of  entropy  on  approach  to  the  eyewall 
((A0e)  =  0.56  K  instead  of  —0.89  K)  with  less  variance 
[cr( A  0e)  =  3.86  K  instead  of  4.72  K],  but  a  cooler  asym¬ 
metric  entropy  at  the  seed  point  ((0')  =  —4.34  K  instead 
of  —3.47  K),  and  this  set  of  trajectories  approaches  the 
eyewall  more  abruptly  (change  in  azimuth  (AA)  =  4.15 
rad  versus  8.06  rad).  Following  the  approach  of  (2) 
above,  the  axisymmetric  mean  eyewall  can  be  found  at 
this  time  to  cool 

A0e(eyewall)  ~  i  p  -  0  27  (+0-56  “  4-34)  K 

«  -  1.40  K.  (3) 

The  degree  of  eyewall  cooling  was  found  to  be  some¬ 
what  more  consistent  with  the  greater  number  of  class 
IV  trajectories. 

4.  Summary  and  conclusions 

This  study  builds  on  the  numerical  study  of  Hurri¬ 
cane  Bonnie  by  Braun  et  al.  (2006)  by  examining  the 
interactions  between  the  various  air  masses  (eye,  eye- 
wall,  inflow,  and  environment)  of  the  TC.  Through 
analysis  of  parcel  trajectories,  we  have  been  able  to 
demonstrate  that  the  low-level  eye  air  can  be  carried 
into  the  eyewall;  that  this  low-level  eye  air  is  replen¬ 
ished  by  inflow  under  the  eyewall;  that  the  mid-  to 


upper-level  eye  air  can  be  exchanged  with  eyewall  air  in 
a  two-way  process;  and  that  mid-  to  upper-level  air  ex¬ 
terior  to  the  eyewall  can  ventilate  the  eyewall  (in  the 
sense  of  Simpson  and  Riehl  1958;  Riehl  and  Malkus 
1961;  Frank  and  Ritchie  2001;  Emanuel  et  al.  2004). 

The  trajectories  examined  demonstrate  that  the  su¬ 
perintensity  mechanism  as  described  by  PM03  and 
Braun  (2002)  is  active  in  the  Bonnie  simulation.  The 
results  suggest  that  the  superintensity  mechanism  is  still 
able  to  benefit  the  intensity  of  the  storm  while  it  is 
being  sheared  by  the  environment.  Trajectories  origi¬ 
nating  in  the  low-level  eye  and  possessing  relatively 
high  0e  are  stirred  into  the  base  of  eyewall  updrafts  and 
introduce  an  additional  source  of  heat  to  the  eyewall. 
The  low-level  eye  is  replenished  by  boundary  layer  par¬ 
cels  that  slip  underneath  the  eyewall  updrafts  to  linger 
for  some  time  in  the  eye  (typically  40-60  min),  thereby 
acquiring  an  elevated  0e  through  prolonged  interaction 
with  the  ocean  surface.  The  stirring  at  low  levels  be¬ 
tween  the  eye  and  eyewall  occurs  predominately  left  of 
the  mean  deep-layer  environmental  shear  vector.  Our 
results  are  consistent  with  the  analysis  of  updrafts  in 
this  MM5  simulation  conducted  by  Braun  et  al.  (2006), 
and  are  also  supported  by  the  observational  studies  of 
Molinari  et  al.  (2004)  and  Corbosiero  (2004),  the  model 
results  of  Frank  and  Ritchie  (1999,  2001),  and  the  theo¬ 
retical  analysis  by  Reasor  et  al.  (2004).  The  results  in 
Bender  (1997)  also  exhibited  asymmetries  forced  by 
imposed  vertical  shear,  though  the  rainfall  in  his  study 
occurred  preferentially  upshear  rather  than  left  of  the 
shear  vector. 

More  than  half  of  trajectories  originating  inside  the 
eye  in  the  middle  and  upper  troposphere  are  not  being 
contained  inside  the  eye  in  our  5-h  trajectory  calcula¬ 
tion  period.  More  eye-eyewall  mixing  occurs  in  the  up¬ 
per  troposphere  than  does  in  the  hurricane  boundary 
layer  and  middle  (z  —  5.6  km)  troposphere.  For  mature 
hurricanes  in  low  vertical  wind  shear  conditions,  W98 
hypothesized  that  eye  containment  of  air  parcels  above 
the  eye  inversion  would  be  quite  strong.  Our  case  is 
more  consistent  with  the  weakening  cases  of  W98  and 
Carsey  and  Willoughby  (2005).  Containment  of  air 
within  the  eye  is  compromised  in  our  case  by  stirring  of 
air  from  the  eyewall  into  the  eye,  but  very  little  air  is 
stirred  into  the  eye  from  the  environment  above  the 
boundary  layer.  Although  a  comprehensive  analysis  of 
mass  replacement  to  the  eye  at  mid-  to  upper  levels  was 
not  carried  out,  there  is  evidence  of  eyewall  air  being 
stirred  inward  into  the  eye  at  these  levels  (e.g..  Fig.  5). 

Ventilation  (Riehl  and  Malkus  1961;  Emanuel  et  al. 
2004)  is  illustrated  in  the  Bonnie  simulation,  as  trajec¬ 
tories  track  the  migration  of  environmental  air  with  low 
6e  into  the  eyewall.  The  net  impact  is  about  a  1-K  re- 
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duction  of  eyewall  average  0e  based  on  analysis  of  air  at 
z  =  5  km  in  the  eyewall  for  two  instances  in  time.  For 
this  simulation  of  a  vertically  sheared  hurricane,  the 
intensity  is  plausibly  much  greater  than  E-PI  (simulated 
55  vs  E-PI  38.2  to  52.2  m  s-1).  The  interaction  of  envi¬ 
ronmental  air  with  low  entropy  and  the  eyewall  is  hy¬ 
pothesized  to  weaken  the  TC.  This  is  supported  by  the 
modeling  study  of  Frank  and  Ritchie  (2001),  in  which 
eddy  mixing  between  the  eye  and  eyewall  was  hypoth¬ 
esized  to  result  in  ventilation  of  the  eye  in  the  upper 
troposphere  and  subsequent  weakening  of  their  mod¬ 
eled  TC. 

One  caveat  to  consider  when  interpreting  the  results 
presented  here  is  that  the  mesoscale  model  used  in  this 
study  is  not  capable  of  resolving  all  the  physical  pro¬ 
cesses  relevant  to  the  eye-eyewall  mixing  problem  in 
TCs.  Subgrid-scale  diffusion  parameterizations  used  in 
most  numerical  models  are  not  applied  with  any  real¬ 
istic  expectation  of  representing  subgrid-scale  turbu¬ 
lence  in  the  atmosphere,  rather  they  are  used  as  a 
means  to  maintain  numerical  stability  in  the  model. 
This  particular  simulation  using  MM5,  for  example,  em¬ 
ploys  the  Blackadar  boundary  layer  scheme,  which  pa¬ 
rameterizes  subgrid-scale  turbulence  (in  addition  to 
boundary  layer  mixing)  on  the  vertical  structure  only 
(MM5  does  not  include  any  horizontal  mixing  terms  in 
the  boundary  layer  scheme;  lateral  mixing  occurs  only 
through  numerical  diffusion).  One  thus  cannot  expect 
the  model  to  realistically  represent  subgrid-scale  mixing 
that  would  occur  in  the  real  atmosphere,  and  this  ca¬ 
veat  should  be  kept  in  mind  when  interpreting  the 
Lagrangian  results  presented  in  this  study.  Addition¬ 
ally,  although  the  2-km  inner  grid  simulation  is  capable 
of  providing  a  representation  of  the  mesoscale  eddies 
occurring  in  the  TC  inner  core,  this  horizontal  resolu¬ 
tion  is  too  coarse  to  represent  the  larger  turbulent  ed¬ 
dies  of  the  Kolmogorov  inertial  range.  Using  a  model 
with  horizontal  grid  spacing  of  500-2000  m  would  still 
not  suffice,  since  the  effective  Reynolds  number  (as 
defined  in  Bryan  et  al.  2003)  of  the  flow  at  this  resolu¬ 
tion  would  not  be  large  enough  to  allow  an  explicit 
calculation  of  the  turbulence.  The  model  configuration 
used  in  this  study  is  the  current  state  of  the  art,  how¬ 
ever,  and  the  authors  look  forward  to  studies  similar  to 
this  one  being  carried  out  on  models  capable  of  resolv¬ 
ing  the  TC  inner-core  mixing  processes.  Computing 
power  has  advanced  to  the  point  where  a  large  eddy 
simulation,  with  a  horizontal  grid  resolution  of  200  m  or 
finer,  could  realistically  be  implemented  in  a  TC  simu¬ 
lation  such  as  the  one  presented  in  this  study. 

Future  work  will  consider  several  of  the  questions 
raised  herein  in  more  detail.  Specifically,  how  substan¬ 
tial  is  1-K  cooling  in  the  eyewall  in  reducing  the  inten¬ 


sity  of  the  TC?  A  continued  study  of  sheared  storms  in 
an  idealized  setting  is  ongoing  to  examine  the  efficiency 
of  ventilation  for  various  magnitudes  of  shear.  Finally, 
it  remains  to  identify  the  principal  asymmetric  struc¬ 
tures  around  the  storms  that  contribute  the  most  to  its 
ventilation. 
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APPENDIX  A 

Determination  of  Storm  Center 

The  technique  for  finding  the  storm  center  described 
below  basically  attempts  to  place  the  maximum  radial 
gradient  in  surface  pressure  close  to  the  same  radius 
from  the  center  of  the  storm  at  all  azimuths.  The  alter¬ 
native  method  of  using  the  location  of  minimum  surface 
pressure  was  found  to  track  centers  of  mesovortices 
rather  than  the  storm-scale  center  of  pressure.  The 
technique  described  below  is  designed  to  avoid  any  ar¬ 
bitrary  selection  of  scale  (the  implicit  scale  is  the  spatial 
scale  over  which  the  signature  of  the  gradient  of  surface 
pressure  is  distributed)  and  any  iterative  processes  for 
solution. 

The  method  for  finding  the  storm  center  is  as  follows. 
Consider  each  grid  point  of  a  select  subdomain  of  an 
x-y  grid  of  simulated  surface  pressure.  Consider  each  of 
four  rays  radiating  along  the  four  cardinal  directions 
from  each  point.  Obviously,  grid  points  near  the  edge  of 
the  defined  grid  will  lack  data  in  one  or  two  directions. 
For  this  reason,  consider  our  select  subdomain  to  be  the 
middle  l/9th  of  the  domain;  that  is, 


2  Nx 

~Y 

(Al) 

2  Nv 

> 

3  ’ 

(A2) 

where  Ix  and  Iy  are  the  index  number  of  a  grid  point  in 
the  x  and  y  directions,  and  Nx  and  N y  are  the  total 
number  of  grid  points  in  the  x  and  y  directions,  respec¬ 
tively.  Determine  the  minimum  M  of  the  number  of 
grid  points  lying  along  each  ray.  Along  each  ray,  a  pro¬ 
file  of  surface  pressure  P;  can  be  found,  where  j  e  {1,2, 
3,  4}  is  an  index  for  each  of  the  four  directions.  Deter¬ 
mine  a  score  for  each  grid  point  Sxy  by  a  least  mean 
squares  difference  from  their  mutual  average 
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M 

EE  Pj(k)  ~  P(k) 

j  k=  0 


where 

S*, 

P  =  (A4) 

is  the  averaged  pressure  profile.  Select  the  point  where 
s  is  a  minimum  as  the  storm  center. 

This  method  works  because  we  know  the  innermost 
grid  of  the  numerical  model  tracks  the  movement  of  the 
central  surface  pressure  minimum,  thus  we  are  assured 
that  the  center  will  be  found  within  the  middle  Vs>  of  the 
domain.  Also,  the  scale  of  the  vortex  is  somewhat  larger 
than  the  scale  of  the  inner-mesh  grid.  When  the  vortex 
scale  is  much  smaller  (such  as  on  the  coarsest  grid  mesh 
of  our  model  run)  than  the  size  of  the  grid,  then  the  four 
profiles  might  easily  miss  the  vortex  entirely  and  the 
score  will  measure  differences  in  the  more-or-less  flat 
environmental  surface  pressure  field,  providing  a  small 
score  trivially. 

We  may  take  M  above  to  be  a  fixed  constant  suffi¬ 
ciently  small  to  handle  every  possible  profile;  that  is, 

(Nx  N  \ 

M  =  mint  —  ,  — j,  (A5) 

but  we  found  that  change  did  not  alter  the  resulting 
centerfind  for  our  case. 

APPENDIX  B 

Computation  of  Trajectories 

A  trajectory  is  the  computation  of  the  path  of  a  hy¬ 
pothetical  fluid  particle  through  a  specified  flow  field 


enough  to  contain  the  complete  translation  of  the  grid 
as  it  follows  the  storm  is  created,  and  all  information  is 
remapped  to  this  grid.  Naturally,  much  of  this  universal 
grid  is  empty,  and  no  effort  is  made  to  populate  these 
empty  areas  with  other  information  available  to  us. 
Any  trajectory  leaving  the  inner  (moving)  mesh  of  the 
MM5  simulation  is  simply  identified  and  considered  to 
have  left  the  near-storm  environment. 

The  grid  stagger  of  the  wind  fields  is  eliminated  by 
interpolating  the  u,  v,  and  w  wind  fields  to  the  0  or 
thermodynamic  grid  points,  which  is  where  the  univer¬ 
sal  grid  is  defined.  The  winds  are  translated  from  m  s-1 
to  the  units  of  (grid  points)  s_1,  taking  into  account  the 
effects  of  latitude  using  the  formulas 

u  u 

“grid  =  At?  =  cos(<(>)|  AA|re  ('B2) 


V  V 

”grid  =  Ay  =  (A<j>)re  ’ 


(B3) 


where  «grid  and  ugrid  are  wind  speeds  in  grid  units,  <f>  and 
A  and  latitude  and  longitude,  respectively,  converted  to 
radians  and  r(l,  is  the  mean  radius  of  the  earth.  Calcu¬ 
lation  in  terms  of  grid  spacing  per  unit  time  is  a  coding 
convenience  that  allows  a  general-purpose,  generic, 
well-tested  trajectory  routine  to  be  used  without  being 
required  to  manage  details  of  spherical  geometry.  If  the 
geographic  pole  were  included,  we  can  expect  the  tra¬ 
jectory  code  to  fail  there  unless  a  more  appropriate 
model  grid  system  were  used. 

Since  the  trajectories  are  computed  in  earth-relative 
coordinates,  storm  motion  is  not  subtracted  from  the 
winds,  but  the  resulting  trajectories  can  be  renavigated 
to  storm-centered  coordinates.  This  renavigation  was 
carried  out  for  all  trajectories  shown  in  this  study. 
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dt 


V(r„„  t). 


(Bl) 
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